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Abstract 

We present a Newtonian multi-fluid formalism for superfluid neutron star cores, focussing on the 
additional dissipative terms that arise when one takes into account the individual dynamical degrees 
of freedom associated with the coupled "fluids". The problem is of direct astrophysical interest 
as the nature of the dissipative terms can have significant impact on the damping of the various 
oscillation modes of the star and the associated gravitational-wave signatures. A particularly 
interesting application concerns the gravitational-wave driven instability of f- and r-modes. We 
apply the developed formalism to two specific three-fluid systems: (i) a hyperon core in which both 
A and hyperons are present, and (ii) a core of deconfined quarks in the colour-flavour-locked 
phase in which a population of neutral kaons is present. The formalism is, however, general and 
can be applied to other problems in neutron-star dynamics (such as the effect of thermal excitations 
close to the superfluid transition temperature) as well as laboratory multi- fluid systems. 
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I. INTRODUCTION 



Neutron stars provide unique laboratories for the study of the state of matter under 
extreme conditions. Comprising roughly one and a half solar masses within a radius of 
about 10 km, these very compact objects are likely to have core densities reaching several 
times the nuclear saturation density. Thus, they offer the opportunity to probe the cold, 
high density region of the QCD phase diagram (that cannot be explored in laboratory 
experiments) . Improvements of our theoretical description of these systems are crucial. Not 
only is a wealth of data from X-ray and radio observations available already, but ground- 
based gravitational- wave detectors have reached their initial design sensitivities and are now 
being upgraded to the second generation level [1] . This is significant since gravitational- wave 
observations would provide truly complementary information about these exotic objects [2]. 
Gravitational waves carry an imprint of the dynamics of the internal, high density, regions 
of the star while the electromagnetic signature relies on a complex interaction with, and 
processes within, the star's magnetosphere. 

Accurate modelling of neutron star interiors is important, but fraught with difficulty. A 
detailed model of neutron star dynamics must not only account for exotic states of mat- 
ter at high density. It must also consider the interaction of the crustal nuclei (the outer 
kilometer or so) with the star's core, the presence of a strong magnetic field and various 
superfiuid/superconducting states. Superfluidity, in particular, adds dimensions to the prob- 
lem by introducing new dynamical degrees of freedom. This can have profound consequences 
for the dissipation in the interior of the star [3, 4]. In the last few years, the effect of super- 
fluidity on the damping of the gravitational-wave driven instability of the so-called r-modes 
has sparked much interest (see [5-7] for recent discussion), and we now know that the details 
of the interior microphysics can play a significant role. Effects such as the vortex-mediated 
mutual friction damping [7-9] and the enhanced bulk viscosity due to the presence of hyper- 
ons [7, 10-12] or of a quark condensate in the colour-flavour locked (CFL) phase [13] can all 
modify the damping of the unstable modes considerably. It has, in fact, recently been shown 
that measured neutron star spins and inferred temperatures in Low-mass X-ray binaries are 
not consistent with the standard models for the r-mode instability window. These results 
suggest that some form of enhanced damping is required [14, 15]. Further theoretical inves- 
tigation is clearly necessary if we want to be able to make direct, quantitative, comparisons 
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with the X-ray data. 

The aim of this paper is to construct a general framework for dissipative multi-fluid 
neutron star cores, building on the variational formulation of two-fluid neutron star hydro- 
dynamics [4, 16] (see also [17, 18]). We will focus on the dynamics of three-fluid systems at 
finite temperature. This problem is relevant both in the case of hyperon cores, where one 
expects not only neutrons and protons, but also a population of A and E~ hyperons [19] , 
and in the case of cores of deconfined quarks in the CFL phase, where one should account 
not only for the quark condensate but also for the thermal excitations of the system and 
possibly a population of kaons [20-22] . Moreover, given the strong density dependence of 
the various superfluid pairing gaps [23] there will always be regions in the star that are 
close to the critical superfluid transition temperature. In these regions thermal effects are 
not negligible. Again, this situation can be described by a three-fluid system comprised of 
neutrons, protons and thermal excitations [24] . 

We pay particular attention to the dissipative terms in the hydrodynamical equations. It 
is well known from the study of superfluid ^He that there are more dissipation coefficients in 
a superfluid system than in the standard Navier-Stokes description. However, although this 
was first pointed out in [3] in the context of neutron star physics, and several authors have 
studied the effect of superfluidity on the standard dissipation coefficients (see [25, 26] for 
recent examples), very little effort has been made to understand the nature (and dynamical 
role) of these extra coefficients. Having said that, there have been some recent attempts to 
quantify the effect of superfluidity on dissipative neutron star oscillations [6, 27]. In this 
paper we build on the work of [4], correcting some conceptual issues, in order to develop a 
general formalism for dissipation in multi-fluid systems. Wc arc developing the formalism 
with neutron star problems in mind, but the framework is sufficiently general that it can be 
applied to entirely different systems, e.g. ones that can be studied in the laboratory. In fact, 
a three-fluid model is required to account for the interactions between phonons and rotons 
(and the associated thermal conductivity) in ^He [28, 29]. A recent discussion of extended 
thermodynamics and causal heat conduction [30] in relativistic models should also be noted. 
Finally, it is worth pointing out the close connection between the multi-fluid model that we 
advocate and the general framework of extended irreversible thermodynamics [31]. 

The paper is structured as follows: In section II we discuss the key length-scales associated 
with, and the general notion of, a "multi-fluid" system. In section HI we outline our flux- 
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conservative multi-fluid formalism, leading to the general form for the dissipation coefficients 
presented in section IV. In section V we discuss the irrotational constraint for superfluid 
flows. Section VI then deals with perturbations of the multi-fluid system for backgrounds 
representing slowly rotating stars, while section Vll is devoted to two explicit examples, (i) 
that of a hyperon core and (ii) that of a core of deconfined quarks in the CFL-K*^ phase. 
Finally, we present our conclusions in section VIII. 

Throughout the paper we use a coordinate basis to represent tensorial relations. We thus 
distinguish between contravariant (of the form v^) and covariant (of the form Vi) objects, 
and raise and lower indices with the (three dimensional flat-space) metric gij. We shall also 
identify the different fluids with constituent indices, so that t'^ will be the velocity of the 
x-th fluid, while relates to the y-th. We do not assume that the constituent indices have 
any geometric meaning (although it is worth noting that it is possible to develop the idea 
of "chemical covariance" ) , so they can be placed either as sub- or superscripts, depending 
on which makes the relevant expression the clearest. 

II. THE DIFFERENT SCALES IN THE PROBLEM 

Before we discuss specific models, it is useful to consider the "big picture" of fiuid dy- 
namics for ultra-dense matter. This is a key issue, because many neutron star scenarios 
require an understanding of dynamics on the macroscopic scale, for which hydrodynamics 
is the natural tool. At the same time, these models must build on an understanding of 
the microsphysics, for which we need to turn to nuclear and particle physics. A properly 
formulated theory of fluid dynamics tracks the evolution of fluid elements — their trajec- 
tories through spacctimc and how the various thermodynamic properties (i.e. an equation 
of state with relatively few parameters) change. Each fluid element is small enough that 
it can be considered as a "point", but large enough to contain many particles so that a 
"smooth-averaged" thermodynamic description is appropriate. 

Let us consider some of the issues involved; first of all, nuclear and particle physics models 
tend to assume global Lorentz invariance for all (matter) fields. This means that gravity 
is neglected. However, one thing that general relativity does well is break global Lorentz 
invariance, while maintaining it locally. Glendenning [32] points out that over the scale of a 
fiuid element in a neutron star, the change in the spacetime metric is negligible and one can 
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always erect a local reference frame where Lorentz invariance holds. Basically, there is a clear 
separation of scales and one can use Lorentz-invariant nuclear and particle physics models 
provided one accepts that "global" in this context means "on the scale of a fluid element" . 

In essence, for a compact star, wc do not have a global phase space to be applied to all 
particles; rather, each fluid element has its own phase space that applies to the particles 
within that element. 

This has repercussions when we consider the micro-physics. For example, in the case of 
fermions (like neutrons and protons) one may ask to what extent one can ignore particle 
states "below" the Fermi surface. On the one hand, one might argue "yes" because the 
lowest energy is zero, and particles occupying that and the other energy states below the 
Fermi surface are locked in, on average, because nearby states will be occupied, and the 
available energy reservoir may not be sufficient to launch particles above the Fermi surface. 
In this line of reasoning, it will only be states near the Fermi surface that participate in 
the fluid dynamics. To some extent this argument is correct, and only particles near the 
Fermi surface contribute to the transport coefficients required in the fluid model (viscosities, 
thermal conductivity etcetera). However, at the same time all particles contribute to the 
global dynamics, as represented by the star's large scale oscillation modes. This is easy 
to see if we consider the fact that these modes couple to the gravitational field, which is 
sourced by all the matter. The fluid dynamics model requires information about both bulk 
properties and transport phenomena, making its formulation a challenge. 

Next, let us consider the scales associated with fluid dynamics. This problem is central 
to the analysis in this paper. In order to discuss a "multi-fluid" system, we obviously need 
to have some understanding of what this concept entails. For ordinary matter, the scale 
is simply set by interparticle collisions. Since we need to associate a single "velocity" with 
each fluid element, the particles must be able to equilibrate in a meaningful sense (e.g. have 
a velocity distribution with a well deflned peak) . The relevant length-scale is the mean- free 
path. This concept is closely related to the shear viscosity of matter. In the case of neutrons 
(which dominate the outer core of a typical neutron star) we would have 



where vp is the relevant Fermi velocity and we have used the estimate for the neutron- 
neutron scattering shear viscosity rj from [25]. 




(1) 
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This estimate gives us an idea of the smallest scale on which it makes sense to consider 
the system as a fluid (about a micron). It also hints at systems with distinguishable multi- 
fluid behaviour. Consider a system with two particle species, and assume that the mean-free 
path associated with scattering of particles of the same kind is (for some reason) significantly 
shorter than the scale for cross-species collisions. Then we have two clearly defined "fiuids" 
and it makes sense to consider the problem using the machinery that we will discuss later. 
This is, however, not quite the situation that is motivating the present work. Our focus 
is on systems that exhibit superfluidity. At the most basic level, superfluidity implies that 
there is no friction impeding the flow. Technically, the mean-free path diverges and the 
previous argument does not work anymore. However, a superfiuid system has a different 
scale associated with it; the so-called coherence length. The coherence length arises from 
the fact that a superfiuid is a "macroscopic" quantum state, the fiow of which depends on 
the gradient of the phase of the wave-function (the so-called order parameter). On some 
small scale, the superfluidity breaks down due to quantum fluctations. This scale is known 
as the coherence length. It can be taken as the typical "size" of a Cooper pair in a Fermionic 
system. On any larger scale the system exhibits fluid behaviour (in the sense of the Landau 
two-fluid model for Helium [28, 29]). For neutron-star superfiuids, the coherence length is 
of the order of tens of Fermi [33, 34]; much much smaller than the mean- free path in the 
normal fluid case. This means that superfluids can exhibit extremely small scale dynamics. 
Since a superfluid is inviscid, superfluid neutrons and superconducting protons (say) do not 
scatter (at least not as long as thermal excitations can be ignored) and hence the outer core 
of a neutron star requires a multi-fluid treatment. 

It would seem then, that one can meaningfully take the fluid elements to have a size 
of the order of the coherence length. However, in reality yet another length-scale needs 
to be considered. On scales larger than the Debye screening length, the electrons will be 
electromagnetically locked to the protons, forming a charge-neutral conglomerate that does 
exhibit friction (due to electron-electron scattering). Moreover, at finite temperatures we 
need to consider thermal excitations for both neutrons and protons, making the problem 
rather complex. 

Furthermore, ideal superfiuids are irrotational and neutron stars are not. In order to 
mimic bulk rotation the neutron superfiuid must form a dense array of vortices (breaking 
the superfiuidity locally). This brings another length-scale into the picture. In order to 
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develop a useful fluid model, we need to average over the vortices as well. This makes the 
effective fluid elements much larger. The typical vortex spacing in a neutron star is of the 
order [33] 

/ p \i/2 

(in « 4 X 10"^ cm , (2) 

\ 1 ms / 

where P is the neutron star spin period. For a slowly rotating object the fluid elements we 
consider (at the end of the day) may be quite large (although obviously still much smaller 
than the size of the star). 



III. FLUX-CONSERVATIVE MULTI-FLUID MODEL 



We take as our starting point the flux-conservative multi- fluid formalism developed in [4] . 
This model combines the conservation laws for mass, energy and angular momentum with 
the variational approach developed by Prix [16] (see also [17, 18]). In this model (which 
represents the Newtonian limit of Carter's convective variational formalism in relativity 

[35, 36]) the main variables are the particle fluxes = n^vf (where is the particle 
number density of the x-th component, and vf its velocity) and the equations of motion are 
derived from a Lagrangian density C of the form 



2n. 

X 

This allows us to define the conjugate momenta 



Px = 



dnl 



08 



dnl 



(3) 



(4) 



where we still need to provide the energy density functional £, which includes the internal 
energy. Following [4] we consider an energy functional that is manifestly isotropic and 
Galileian invariant, i.e. we take 

S^S{n^,wl^), (5) 

with the relative velocities defined by 



w. 



VI - v:, 



and 



The chemical potentials are then obtained from 

d£ 
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(6) 



(7) 



while the entrainment coefficients are given by 

88 



(8) 

nx 

These coefficients describe the fact that the momentum of each species is not necessarily 
parallel to the associated flux. Instead, it takes the form: 

Trf = tv^pI = m^n^vf + 2 ^ "yx^y^ ■ (9) 

The most common context in which this effect has been studied relates to neutrons and 
protons in neutron star cores. In this case the effect is due to the fact that, because of 
the strong nuclear interaction, each neutron is associated with a virtual cloud of protons. 
This modifies the effective neutron mass in a dynamical setting [16, 37]. Recently, the 
entrainment parameters in a hyperon core have also been discussed [38, 39] . The usefulness 
of the entrainment concept for entropy, and its relation with thermal relaxation and causal 
heat conduction has also been explored in recent work [30, 40-42] 

The entrainment is a dynamical effect that arises naturally within the variational model. 
In a practical application, it depends on the nature of the microphysics that one includes in 
the equation of state (represented by £). As discussed in [4] the momentum equations take 
the form: 

dtK + Vj (4< + A'O + nxV, (^/.. - ^m^vl^ + p.V,$ + TTjV.vi = ft , (10) 

where the tensor Df^ represents the viscous stresses, while the "forces" ff allow for the 
transfer of momentum between the two components. The particle mass densities are px and 
$ is the gravitational potential. The latter satisfies the usual Poisson equation 

V'$ = 47rG'J]px. (11) 

X 

In the following we shall consider the case of an isolated system, for which ff — Q. 
The continuity equations can be written 

dtn^ + V,- {n^vi) = T^. (12) 

We will assume overall mass conservation for our system, which in terms of the particle 
creation rates leads to 

5^mxrx = 0. (13) 



This constraint makes sense in Newtonian theory, but is not necessarily justified from a 
microphysics point of view, as it is baryon number, rather than mass that is conserved in a 
given nuclear reaction. 

Finally, the total energy creation rate per unit volume is given by 

1 



^ext 



</r + Dt,V^< + Ux - 



(14) 



So far the model is quite general, but we now want to make direct contact with thermo- 
dynamics. In doing so, it makes sense to highlight the entropy component. We will, after 
all demand that the various dissipation channels adhere to the second law. For an isolated 
system we have e^^* = 0, so that if we consider one of our fluids to represent the entropy of 
the system [66] (we denote this component by the subscript s), the above relation leads to 

TV, = -f,vl - D,,VX - E (^x)"^ + /r< + V^<) . (15) 

where 

fi = X^x/i" total force acting on the system, and Dij = "^^Dfj. Note that, as wc are 

considering a closed system in the following, we take fi = 0. Not also that, the individual 
Dfj do not have to be symmetric, even though the sum D^j must be. Finally, the entropy 
component is taken to be massless, and its chemical potential is the temperature /is — T. 



IV. THE ONSAGER APPROACH 



Let us now move on to consider the general form of the dissipation coefficients. To do this 
we continue to follow the analysis in [4] , and make use of the Onsager symmetry principle 
[43]. Our analysis will, however, differ from that in [4], as in that work the authors neglected 
(erroneously) a number of terms involving Vji'f . Hence, the resulting dissipative model was 
not as general as it should have been. This has already been noted in [28], where the analysis 
was reworked in the case of two fiuids, one of which represented the entropy. As we intend 
to apply the discussion to the case of three (or more) fluids, we flrst consider the general 
form of the dissipation coefficients, turning to particular examples later. 
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For any system, perturbations of the entropy density s away from equilibrium must be 
given by quadratic deviations. This allows us to write [4] 

a,b 

which can be expressed, in terms of the entropy creation rate Fg, as 

1 ^ 
TF, = -- J] XaL-'X, = J2 J"Xa , (18) 

a,b a=l 

where the are known as "thermodynamic forces". They represent a measure of the 
departure of the system from equilibrium, while the "thermodynamic fluxes" 

J« = -i^L«%, (19) 

b 

arise in response. The Onsager symmetry principle states that microscopic reversibility 
implies L"^ = L^". 

By comparing equation (18) to equation (15) we can, by constructing the most general 
form for the tensor L"'' in terms of the thermodynamical forces in the system, obtain the 
most general description of the dissipative terms in the Euler equations; the stress tensor 
Dfj, the forces /f and the terms associated with the reaction rates. With this in mind, let 
us write down the most general form for the tensor L"''. 

From equation (15) it might be tempting to take the thermodynamic forces to be u^xs, 
VjW*g, Vjvl and /Xx, as in [4]. This is not quite appropriate, however, as we need the forces 
to vanish when thermodynamic equilibrium is reached. Hence, we should not work with 
the chemical potentials, which obviously do not vanish in equilibrium. This point comes 
to the fore when we consider problems with reactions, as in the case of bulk viscosity in 
a multi- fluid setting [6]. We need to replace the chemical potential with a more suitable 
"force" , and the affinity [44] is the natural choice. In the context of neutron stars, this point 
has been made by Carter and Chamel [45]. The following analysis combines the key points 
of their discussion with the general multi-fluid formalism from the previous section. 

Suppose there are N total reactions among M various constituents x of our multi-fluid 
system, to be characterized in the usual way as stoichiometric relations between the particle 
number densities [67] u"^ — n^/ (J^x^^) ' 

M M 

^R^^^x^^p/^x ^ / = i,...,iv, (20) 

X X 
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where and are, respectively, the reactant and product stoichiometric coefficients. The 
affinity of the 7*^ reaction is then defined as 

M 



At thermodynamic equihbrium the affinities vanish, which is why they make appropriate 
thermodynamic forces. 

While it is clear that the affinities provide a natural description of the problem, it is 
important to recognize that the formulation is not quite complete at this point. In partic- 
ular, the chemical potentials /i^ become somewhat ambiguous in a multi-fluid context. In 
principle, the chemical potentials should be defined as the energy per particle in the ref- 
erence frame where the chemical (or nuclear) reactions occur [18]. However, a multi- fluid 
mixture is characterized by the presence of distinct velocity flelds, neither of which provides 
the required frame. The relevant frame may, in fact, not be known a priori. The formula- 
tion we consider assumes an expansion away from "equilibrium" , which ultimately involves 
both dynamical and chemical considerations. The equilibrium frame may well depend on 
the dynamical evolution of the whole system. This complicates the issue considerably. Of 
course, in many situations of practical interest this problem may not be too severe. This is 
particularly the case when the relative velocities between the different fluid frames are small 
enough that it makes sense to linearise the problem. Noting that a satisfactory solution to 
the conceptual problem still remains to be developed, we proceed on the assumption that a 
low-velocity model makes sense. 

According to Hess's Law, for each chemical reaction there is only one thermodynamic 
variable to track in order to determine the changes; namely, the "degree of advancement" ^7 
for the various reactants. For each of the / = l...A^ reactions, a variation A^/ corresponds 
to a variation Ai/f of the participating fluids: 




(21) 



X 




'r Az/? Auf AuJ 

- = •••= rT = -^r = - = -^r = ' (22) 

^^S ^ U ^ V 

represent the x-components for which the and P^ are non-zero. 



The (irreversible) change As in the entropy due to these reactions is given by 




(23) 



7=1 
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By comparing with equation (17), we see that the A^/ represent the appropriate thermody- 
namic "fluxes". 

The variations of the individual number densities, in some time interval At, can also 
be determined by 

Ai/^ = r^At , (24) 

where Fx is the particle number creation rate. 

Each of the N reactions has a corresponding change z/j that contributes to Ai/^, with the 
net result (as At — >• 0) 



Hence, 



d 

rx = E(Rx-P3f ■ (26) 

If we define the reaction "velocity" = ^ to be the thermodynamical flux, then the 

change in entropy due to the reactions is 



In the general framework the thermodynamic forces will then be , w^g, VjV*, Vj-w^g, 
and the corresponding fluxes are —Vi,—ff, —D^j and —D^j- Given this, we can follow the 
strategy of [4] to construct the fluxes out of the forces, limiting ourselves to the inclusion 
of quadratic terms. By making use of the Onsager symmetry principle, and noting that 
conservation of angular momentum [68] requires D^j to be symmetric (see equation (22) of 
[4]), we then arrive at; 



-ft 



x,y^s J 

E + 



E 
E 



(28) 
(29) 

(30) 

(31) 
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with 



-^x — 7x — 7x ^ij — ^ 9i3 ^ij — ^ 9ij J 



(32) 
(33) 



Lljki = C9ij9ki + (9ik9ji + 9ii9jk - '^9ij9ki^ , 
^jki = C^9ij9ki + V""^ [9ik9ji + 9ii9jk - \gijgki] + \(^''^eiimt"'ki 



= C^^giogki + 9ik9ji + 9ii9jk - o9ij9ki + 7;(^^''^ijme"^ki ■ 





(36) 



(34) 



(35) 



These relations suggest that the complete set of dissipation coefficients is given by 



essentially the same as in [4]. The important difference is that, not only have we now 
correctly identified the affinities as the thermodynamic forces, but we also have a host of 
new terms relating to the gradients of the entropy velocity vl- These were neglected in [4]. 
An already complex problem has been made richer... 

The problem may seem almost unmanageable at this point. In particular, how are we 
going the determine all the different dissipation coefficients? The answer is that we need to 
complement the phenomenological multi-fiuid model with insights from microphysics (e.g. 
kinetic theory). Once we turn to that question we see that the complexity of the problem 
may reduce drastically. In fact, it is easy to argue that many of the different coefficients in 
the model will have the same microphysics origin (although various geometrical factors may 
differ). To see this, let us consider an example. Let us ask how particle scattering enters the 
problem. In a single-component system, the answer is simple. Scattering leads to friction 
that prevents fluid clement from shearing, i.e. determine the coefficient of shear viscosity. 
The multi-fluid setting is more complicated, because we have to account for scattering both 
between particles of the same species and inter-species scattering. The former will (again) 
lead to the familiar shear viscosity from the Navier-Stokes equations, one term for each 
fluid species. The inter-species scattering affects the relative degrees of freedom. Two 
fluids can flow linearly through one another, and they can also have relative shear and 
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expansion. The corresponding coefficients will all relate back to the scattering rates. A 
similar analysis relates to the various bulk viscosities. In this case, the problem reduces 
because the dissipation originates either from the relevant reaction rates (the case of main 
importance for neutron stars) or from fluctuations of the internal degrees of freedom for each 
species (as in water). A useful example, with direct relevance for one of the models discussed 
below, has been analysed by Gusakov and Kantor [10] (although it should be noted that a 
translation between their model and that considered here is non-trivial). Other interesting 
discussions can be found in [46] and [47]. 

One may also simplify the problem by constraining the physics. In the present con- 
text, the most relevant constraint is associated with superfluidity, and we now turn to the 
corresponding problem. 



V. THE SUPERFLUID CONSTRAINT 



Let us now assume that one of the fluids (labelled S in the following, not to be confused 
with the entropy component which is represented by a lowercase s) is superfluid. In this case 
we would expect it to be, at least in the limit of low temperatures and velocities, irrotational. 
Following [16] and generalising the work of [28], we impose the irrotational constraint on 
the momentum of this fluid. This means that we require 



e'^'V^pl = , 



(37) 



which leads to the constraint 



V.^ = -iff - VW%) , 



ns 



(38) 



for some scalar In order to analyse this constraint it is useful to rewrite equation (15) in 



terms of the variables jx 



nxW*g. This leads to 



where 



TT-x 



1] 



and 



1 



/r- 



The constraint in (38) thus takes the form 



nx 



(39) 



(40) 



(41) 
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Repeating the analysis of the thermodynamic fluxes (from the previous section) in terms of 
the new variables, we flnd that 







(42) 




— ^^^ij -Jys ^ ^^^ijk ^ Jys ' 

= E E ^^^^ + M^JklV'vi + M^,, V% , 


(43) 


-Dij 


(44) 




= E E ^^'^^ + + Mr.^yvi + m->,, v% . 


(45) 



(where the M coefficients essentially replace the L's from the previous section). The con- 
straint in (38) now implies that, for the superfluid component one must have 

V,* = (^Mf/Aj + M%yi^ + Mf^.iS/'^vi) + (y^Mfy, - Mf^) , (46) 

which implies 

Mf^ = ^'Mf^ = • (47) 

To see that these coefficients must vanish individually, consider Eq. (33). The coefficient Mf^ 
has the same form as L^^ in Eq. (33), i.e. it is symmetric, while the Mfjf. is antisymmetric 
when i and k are interchanged. 

The superfluid constraint thus takes the form 

V,* = [Mf/Aj + M%iV% + M!,,yvi) , (48) 

from which we infer that 

M^' = r^9ij , (49) 

= C^'grni , (50) 

Mfjki = C^9ij9ki ■ (51) 

If we compare these results to the general expressions in (36) we see that imposing the 
superfluid constraint on one, or more, of the fluids signiflcantly reduces the number of 
dissipation coefficients. 

It is important to appreciate that this reduction comes about once we identify the ap- 
propriate thermodynamical fluxes in the system. This highhghts the fact that there may be 
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different ways of formulating any given problem, possibly leading to systems of seemingly 
different complexity. This is somewhat unfortunate, but we are not aware of any general 
prescription for avoiding this ambiguity. It is also worth noting that it is generally not 
meaningful to "translate" the coefficient in the two models we have provided. The basic 
reason for this is that the two systems are "separated" by a phase-transition (at the critical 
temperature for superfluidity). The models we have outlined apply either above, or well 
below, the relevant temperature. The increasing role of thermal excitations as the transition 
is approached make the corresponding problem tricky. In principle, the system must "switch 
on" the various dissipation channels that were removed by the superfluid constraint as the 
critical temperature is approached (reverting to the general system discussed in the previous 
section). 

The reduction associated with the superfluid problem may only be of formal interest, 
however. In practice, the irrotational constraint is too severe since a superfluid can rotate 
by forming an array of vortices. This complicates the problem rather than simplifying it. 
First of all, we need to use the expressions in (36). Secondly, we need to worry about 
additional dissipation channels that come into play when vortices are present. Having said 
that, there may be some merit to an argument that one can take the irrotational model as 
starting point, adding only the particular mechanisms that are due to the vortices. In such 
a model, which may turn out to be accurate in many situations, the reduction of complexity 
due to the irrotational constraint is helpful. Most current discussions of neutron star vortex 
dynamics build on this idea. 

The main problem with the irrotational assumption is that it means that we eliminate the 
vortex-mediated mutual friction, which arises from the balance between the Magnus force 
and the dissipative drag forces (such as electron scattering of vortex cores) on the vortices 
[48-50]. This is known to be a key mechanism in neutron star dynamics. We will not 
discuss the mutual friction in detail here, but it is nevertheless worth making a few general 
remarks. In the standard two-fluid system of neutrons and protons, and if one considers 
straight vortices, the friction between a superfluid component x and another component y 
can be taken into account by explicitly including a force of the form: 

/f^ = B' p^n^€ijkK^w'^y + Bpyji^eijke^^'^k^ Kiw"^ , (52) 

where is a vector aligned with the vortices, with magnitude k — /i/2mx [69] (a hat 
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indicates a unit vector) and Uy is the vortex number density per unit surface area. B and B' 
are coefficients tfiat encode tfie strength of the mechanism, but only the former is associated 
with actual dissipation. Given this additional force we see that, in the presence of vortices 
our analysis of the dissipation coefficients in the system is formally incomplete. We could 
make our analysis more general by accounting for the preferred direction associated with 
when we design the dissipation terms. This extended model would obviously allow for 
the standard mutual friction represented by (52), but also for the presence of a (significant) 
number of additional dissipative terms coupling the fiows to the vorticity (see [28] for a 
discussion) . 

Before moving on, it is also worth noting that the standard form for the mutual friction 
force may not be entirely appropriate. In many circumstances the supcrfluid flow is expected 
to be turbulent. This means that the vortices are no longer "straight" but form a complicated 
tangle. In this case one can no longer use the expression in (52). Alternative forms in which 
the force is proportional to the cube of the relative velocity have been proposed. The form 
of this turbulent force and the coefficients involved are not well established in the neutron 
star context, although there have been recent attempts to understand the relevance of the 
effect [51-53]. 

In fact, the nature of mutual friction in a general multi-fluid system may be considerably 
more complex. One may have to account for the interaction between several kinds of vortices 
which can form in the different superffuid or superconducting condensates. For example, [54] 
has recently shown that for strong entrainment, or in the presence of superffuid E~ hyperons, 
the usual picture in which one has rotation-induced neutron vortices, but not proton vortices, 
does not hold. Rather one can have "composite" vortices, strongly reducing the interaction 
between the superffuid and the magnetic-ffeld carrying superconducting components of the 
star. Another possible complication has been discussed in [55]. Clearly, more work is needed 
in order to understand how to include such concepts in our general picture. 

VI. PERTURBATIONS 

The development of the dissipative multi-ffuid formalism is obviously somewhat abstract, 
and we need to consider explicit examples in order to sec how it can be applied. Ultimately, 
we are interested in how dissipation affects the modes of oscillation of a rotating superffuid 
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neutron star. Given this, and the fact that dissipation is associated with deviations from 
equihbrium, let us consider hnear perturbations of the multi-fluid equations of motion, (10) 
and (12). Assuming that we are interested in rotational instabilities, hke that of the Coriohs 
driven r-mode [5-7], we work in the slow rotation approximation to linear order, meaning 
that wc perturb around a spherical background in which all fluids co-rotate and are in 
hydrostatic and chemical equilibrium. These assumptions simplify the background equations 
considerably, as there is no dissipation and no terms involving relative velocities [70]. One 
simply has to consider the standard equations for hydrostatic equilibrium: 

ViP + pVi* = 0, (53) 

where $ is the gravitational potential (as before) and where total density and pressure are 
given by 

P = I] Px , (54) 
Vip = ^ rix Vi//x ■ (55) 

X 

The perturbed Euler equations in the rotating frame can then be written in the form (rep- 
resenting Eulerian variations by S) 

where is the angular velocity of the star and px = 'tt.x'T'x- The perturbed momentum 
Sirf will not in general be parallel to the perturbed velocity of the x component; due to the 
entrainment effect it acquires components along the relative flows. We have 

Snf = Qij ^^njvi + 2 ^ a^'^Sw^y^ , (57) 

where the a^^ are the entrainment coefficients defined in (8), the perturbations of which we 

do not need to consider since we are considering a co-moving background. We also need the 
perturbed continuity equations, which take the form 

9t5px + V,(pxK)=^"rx, (58) 

where we have assumed that the background is in chemical equilibrium. In other words, we 
take the reaction rates Fx to arise at the linear perturbation level. This obviously makes 
sense since the reactions are triggered by deviations from chemical equilibrium. 
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The relations in (56) and (58) have the same structure for each fluid and represent the 
equations of motions for the x coupled degrees of freedom. In many situations, it can be 
an advantage to make use of this "symmetry" . At the same time, it is instructive to show 

that the equations can be combined to regain the usual Navier-Stokes equations. To do 
this, we sum the perturbed continuity equations in (58) and assume mass conservation (i.e. 
^^m^Fx = 0). This leads to 

dt6p + Wj{pSv^) ^0 . (59) 
where we have introduced the velocity, v\ associated with the total mass flux; 

Meanwhile, the sum of the Euler equations in (56) leads to: 

pdtSvi + ViSp + 2peijkn^v'' - ^Wip + pViS<^> = -V^Dij , (61) 

P 

where 5Dij = Y^^SDfj. 

Clearly, equations (59) and (61) only account for one of the N dynamical degrees of 
freedom, and one has to supplement them with the equations of motion for the remaining 

— 1 degrees of freedom, either directly from equations (56) and (58) or with suitable 
combinations of these. In most previous work, focused on the two-fluid case [7, 56-58], it 
has been found advantageous to work with the difference of the Euler equations. This leads 
to an evolution equation for the relative velocity, sourced by the deviation from chemical 
equilibrium. In the general case, where one has more fluids and different kinds of reactions 
it is not so obvious what the best combination to use may be. One may have to consider 
the issue on a case by case basis. 

Before we consider explicit examples, it is worth commenting on one particular problem 
area where the present results may be applied. Problems in neutron star dynamics are 
closely hnked to the effort to detect gravitational waves. A key problem concerns oscillations 
and instabilities of rotating compact stars. In order to consider neutron star models with 
realistic interior composition, one needs a formalism able to determine the timescale on 
which oscillations of superfluid neutron stars are damped out by various dissipative processes. 
There are essentially two approaches to this problem. The first consists of solving the full 
dissipative problem, essentially "integrating" equations (56) and (58). Given the complexity 
of the different dissipative terms this approach is, however, often not viable. The second 
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option is to estimate the various dissipation timescales from energy integrals, see [58] for 
discussion. This approach is based on the assumption, vahd in many physical circumstances, 
that the dissipation is weak and does not significantly alter the nature of the solutions to 
the conservative problem (in which the dissipative terms are absent). This will be the case 
when the damping timescale is significantly longer than the dynamical timescales we are 
interested in, i.e. the oscillation periods considered. In this case one "simply" has to solve 
the continuity equations in (58) and Euler equations in (56) without the dissipative terms Df^ 
and force terms f^. In addition, one requires an energy for the system. The relevant object 
can be obtained by multiplying the Euler equations for each component by p^v^ (where 
the * represents complex conjugation) and adding it to its complex conjugate. Combining 
the contributions leads to a total time derivative of a quantity which we can define as the 
energy, E. The time derivative of this energy dtE follows from the "right-hand side" of the 
dissipative equations of motion, and allows us to estimate the damping timescale 

(62) 



dtE • 

Examples of this kind of analysis are discussed in [59] . 



VII. EXAMPLES: EXOTIC NEUTRON STAR CORES 

So far, our discussion has mainly concerned the general multi-fluid formalism. Once 
we move beyond the singlc-fiuid model the situation clearly becomes very complex. The 
main lesson is that we need to consider (at least in principle) a plethora of new dissipation 
channels. In order to gain better insight into this reality, we will consider two problems with 
direct relevance for neutron star astrophysics. Both problems concern the deep neutron star 
core. In the first case we assume that the core has a sizeable hyperon component while the 
second example concerns a deconfined quark core. The examples are similar in that they 
both require us to consider three coupled "fiuids". Yet, they are sufficiently different to 
illustrate the subtleties of these kinds of problems. 



A. Hyperon cores 



The first of our examples concerns a neutron star with a hyperon core. Several proposed 
equations of state predict the appearance of hyperons at supranuclear densities (see e.g. 
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[11, 19]). The recent measurement of a neutron star mass of 1.97 Mq [60] appears to place 
stringent constraints on equations of state with softening components, hke hyperons, but it is 
important to keep in mind that the presence of hyperons is expected for fundamental physics 
reasons [61] and the models remain incomplete until the many-body interactions are fully 
accounted for. As this is an immensely difficult problem, it is relevant to consider indirect 
evidence for the presence (or, indeed, absence) of hyperons in the core of astrophysical 
neutron stars. The r-mode instability may provide interesting constrains in this respect. 

The A and E" hyperons are predicted to have the lowest thresholds for formation. The 
resultant problem is of great interest, in particular for gravitational-wave physics, as the 
presence of hyperons increases the strength of the bulk viscosity and reduces the part of 
parameter space in which gravitational-wave driven instabilities may operate. The exact 
details of the damping can have observational consequences. Two of us have recently ex- 
amined the effect of the additional damping coefficients on the r-mode instabihty [6] . That 
analysis was based on the simpler case of a E" hyperon core, which can be cast as a two-fluid 
problem. In the following we develop present the formalism for a core in which both A and 
E~ hyperons are present. As we will see, this is essentially a three-fluid problem. 

We start from the equations of motion for a cold four fluid system, formed by neutrons 
(n), protons (p), E~ and A (in the interest of clarity, we do not account for the presence of 
muons, even though they would be easily included in the model). This means that we have 
already assumed that the electrons are locked to the protons on a much shorter timescale 
than the dynamical timescales we are considering [71], and neglect their mass compared to 
that of the other components [34, 49]. For similar reasons, the E~ hyperons can be considered 
locked to the charged component. Hence, we let = fp such that all remaining "fluids" 
are overall charge neutral (the scale considered is assumed to be larger than the relevant 
screening length). We will, however, retain a separate E~ mass fraction in the continuity 
equations. We do this in order to make the analysis of the bulk viscosity damping simpler. 
After all, the expectation [11, 62] is that the bulk viscosity will be dominated by hyperon 
creating processes like 

n + n # p + E- , (63) 
n + p ^ p + A. (64) 

Moving on to the perturbations of the hyperon core, we write the momenta of our three 
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components (n=neutrons and c=protons (p) locked to electrons and E hyperons (E), A = A 
hyperons) as 

< = [m^n^vi - 2{a^^wi^ + a^'^wi,. + , (65) 

< = 9ij {(mprip + msns)^J + 2[(a"P + a^^X - , (66) 
vrf - {mAriAT;! + 2[(a^^ + + - (a^^ + a^^XJ } . (67) 

We can then write the perturbed Euler equations as one equation for the centre-of-mass 
velocity, equation (61), and two equations for relative velocities. As there will be no pertur- 
bations of the entrainment for co-rotating backgrounds, the two difference equations take 
the form: 

1 + ^] V/Pf + (-^ ^] V/Pf + -^V,dDi+ 

1-63 PH+PpJ ' VI- £3 PS + Pp/ ' PS + Pp ' 

+ (i + T^ + -^)sJ7-(t^ (68) 

V l-£3 Pn + PvJ VI -£3 Pi: + PpJ 



-(1+ Pn + PE + Pp \ ^ ^^,n ^ A _ ^ PJ^\ ^^^jA 

V ps + Pp y ' V i-£iPE + Ppy ' 

Ps + Ppl-£i 

1-^ P^]5:Ft+(l+ P'^ + P^ + pA sP^. (69) 

In these expressions the tildes (~) indicate that the variable has been rcscalcd with the 
relevant mass m^, e.g. 5pc — ^l^c/fT^^- We have accounted for the different masses of 
the individual components, but it is worth noting that neglecting the mass difference of 
the hyperons, taking all masses to be equal, is a reasonable approximation if we want to 
determine the damping timescale of the r-mode instability [6]. This obviously also makes 
the problem more tractable. 

We have also defined various combinations of the entrainment parameters 

£1 = ^ + ^ ■ '- , £2 = ^ , (70) 

Pn PS + Pp Pn PS + Pp 

£3 = ^ + ^ , £4 = ^ - . (71) 

Pn Pa Pn Pa 
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It is worth noting that, while the problem would simplify significantly if we were to ignore 

these coefficients we have physics justification for doing so. In fact, based on the familiar 

case of neutron- proton entrainment, there is every reason to expect the mechanism to be 

relevant also for the hyperons. This is, indeed, borne out by the results in [38, 39]. 

From the perturbed continuity equations we derive the usual relation for the total density 

perturbations, equation (59), and three equations for the mass fractions 

1 m r 
dt5xp = --Vj {pXp[xA6w^^^ - {1-Xp- xs)Swip]} - Sv^WjXp H ^ , (72) 

dtSxA = {p^a[(^P + ^^)^'^ip - (1 - ^A)'^^nA]} " Sv^VjXA + '^^^^ , (73) 

dtSxj: = {p^A^A^'^iA - i^-^p- 2^s)(^w^p]} - Sv^WjX-E + ^^^^ , (74) 

where we have introduced the fractions Xy^ = p^/ p. Note that the above equations are not 
independent (otherwise one would have four, not three, fluids). We have considered the 
hyperon and proton fractions separately in order to have the hyperon and proton creation 
rates explicit in the equations. The intention is to make the analysis of bulk viscosity simpler 
when one has to deal with the rates of several reactions, as in (64). One then needs to apply 
the additional condition of charge neutrality, which implies that Up — 71^, + rie where Ue is 
the number density of electrons. However, when the fraction of S~ is significant the electron 
population is depleted and we may be able to use the approximate relation Up ^ ns [6]. 

Let us now examine the nature of the dissipative terms in the equations of motion. We 
shall assume that the neutrons and A hyperons are superfluid, and thus apply the irrotational 
constraint (38) to both these components. Meanwhile, following [4], we shall neglect heat 
conduction and assume that the entropy flows with the charged components, which thus 
represents the "normal" component. Note that we are explicitly taking v^, = Vp, so that 
jss = ips = 0. In this approximation the required thermodynamic fluxes take the form 



= E E ^x-'^J + Mf^W.vi + MffV'jl + M^-^VVL , (75) 

X ,/ 

-A,- = E ^i^J + Mm'^H + M-kN'jlp + m4,V% , (77) 

J 

-V- = E ^i'^^ + ^^1^^% + ^Wki'^'lv + Kki'^'Ap ' (78) 
J 

= E ^5'^-^ + ^kN% + V% + M^^, V% , (79) 
J 

23 



with coefficients 

Mi' = -f'J Mi=r'g,j Mf/ = r^' g,, M^f = T^'g,j, (80) 

Mijki = CgijQki + i^ikQn + guQjk - '^gnQki^ , (81) 

= C9^J9kl M^jki = C''9^m , (82) 

M'^'^ = rgim = C^'^i.-^fe^ = C^'^'grm ■ (83) 

We tlius liave tlie standard sliear viscosity coefficient t]^ that appears in the Navier-Stokes 
equations, but there are now 6 bulk viscosity coefficients and the reaction rates depend 
on the ffows via the t"^ and t^"^ coefficients. It is important to note that the new viscosity 
coefficients in the equations of motion are related to the relative flow. A generic perturbation 
of the system, or indeed a typical oscillation mode, will be represented by the coupled degrees 
of freedom. It is by no means obvious from the outset to what extent the various degrees 
of freedom are excited in a given situation. A useful illustration is provided by the (polar) 
f-mode of a superfluid neutron stars, for which the damping is dominated by the vortex 
mutual friction and hence relies on the relative motion of the two components [58]. The 
damping can be strong because the excitation of the relative motion is considerable. This 
is in clear contrast to the results for the (axial) r- modes [7] , for which the exitation of the 
relative flow is weak (at least for the mode that is the strongest gravitational- wave emitter). 

B. CFL and kaons 

Our second example concerns a core of deconfined quarks in the colour-flavour locked 
(CFL) phase combined with a population of kaons. This problem is of great interest as the 
ground state of matter in neutron star cores has been the object of vigorous investigation (see 
e.g. [20] for a review). While the pure CFL phase (which can be conveniently be described 
by a two-fluid model with a CFL condensate coupled to phonons [13]) is the ground state of 
cold matter at asymptotically high densities, it is thought that a kaon condensate is likely 
to be present at realistic neutron star core densities (leading to the so-called CFL-K'' phase 
[21]). This possibility has recently been considered in connection with the damping of the 
r-mode instability [63], but the multi-fluid aspects of the problem have (so far) been ignored. 

In order to model this situation we need to consider a three-fluid system given by the 
(neutral) CFL condensate, the kaons and the excitations of the system, which we treat as 
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a massless entropy fluid (as in [13, 30]). The analysis then proceeds (essentiaUy) as in the 
previous section. We indicate the CFL condensate with c, the kaons with k and the entropy 
with s. The momenta of the different fluids are 



7rf 



(84) 
(85) 
(86) 



Note that, as we impose that the fluids flow together in the background, the unperturbed 
entropy momentum (86) vanishes. From the Euler equations (10) one then flnds that 



sViT = , 



(87) 



which impUcs that the unperturbed core has uniform temperature (we have made the iden- 
tiflcation Ug = s, Hs = T). The background system is isothermaL 

We can again write the perturbed Euler equations as one equation for the total velocity, 
equation (61), and two equations for the counter- moving velocities: 
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_ Qi + Pc(gk -£cs) ^ ^^ij ^ Qi-Pk(gk-£cs) ^ ^^ij 
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+ ^ + - '-^ 5r.-^- P'^^t " , (89) 
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where S(3 = dUc/mc — ^/^k/'^^^k 
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Since the background is co-moving, perturbations of the entrainment parameters do not 
appear in the equations of motion. In contrast is the previous example, most of the required 
coefficients have not been previously considered in the literature. Hence, we do not know if 

it necessary to retain the entrainment between quarks and kaons. The entropy entrainment, 
on the other hand, is known to be key to the thermal relaxation in the problem [30]. The 
effect may be weak at low temperatures, but there are fundamental arguments for why it 
should be present. 

Prom the perturbed continuity equations we derive the usual relation for the total density 
perturbations, equation (59), a continuity equation for the kaon mass fraction (xk = Pk/p) 

dtSx^ = -Vi [xk(l - Xk)(5<] - 5v'WiX^ + ^ , (92) 
P P 

and a conservation law for entropy 

dt5s + Vi{s5vl) = dt6s + Vis{6v' + Xk5<k - <^<J = . (93) 

Here it is worth noting that, since we are working at the linear perturbation level and the 
background is isothermal (i.e. there is no heat flow in the unperturbed configuration) the 
problem is adiabatic. 

Let us now turn to the dissipative terms in the Euler equations. If we assume that the 
CFL condensate and the kaons are superfluid we can apply the irrotational constraint (38) 
to both these components. Meanwhile, the entropy represents the "normal" fluid. The 
situation is thus very similar to that for hyperons, and the fluxes are explicitly given by; 

-V' ^"^^ + ^i^^'^s + M^/V% + Mi'V% , (94) 

X J 

- J-c ^ _jrk ^ , (95) 

-A,- = Yl ^^i'^J + MmV'vi + M^^-fc, V^jL + M;^.,, V^jL , (96) 
J 

-Vt^ = J2 M!/Aj + MyV^vl + Mt^,,V^t + M^^V'jl , (97) 
J 

-Vl = ^t'Aj + V'^s + M^^^'L + V^jL , (98) 
J 
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with 

M," = 7," Mi^T'g,^ Mt^ = r''9^3 ^ t^^' , (99) 

Mijki = CgaQki + rf [gikQji + 9ii9jk - '^Qijgki^ , (100) 

M^jki = CgijQki Mfjki = C^^gimi > (loi) 

Once again we have the standard shear viscosity coefficient rj^ and 6 bulk viscosity coeffi- 
cients. 

In the above analysis we accounted for a charge-neutral superfiuid K'^ condensate. An- 
other possibility is that, for sufficiently hot neutron stars, there exists a sizable thermal 
population of kaons [22, 64]. This would change the problem considerably, as one would 
expect that (at least as a first approximation) the thermal kaons would be locked to the 
entropy (the thermal excitations), thus adding to the "normal" fluid. Moreover, the thermal 
kaons would be massive. To deal with this situation, we could once again solve a two-fluid 
problem, similar to that of a CFL condensate coupled to phonons [13]. The main difference 
would be that the superfluid condensate would no longer be coupled to a massless fluid, but 
rather to the fluid formed by thermal kaons and phonons. The mass density of such a fluid 
would be temperature dependent. This situation has not been considered previously. It 
could be interesting, both from a conceptual and a practical point-of-view. There may, for 
example, be interesting consequences for the r-mode instability, as the Coriolis force which 
drives the mode will not act on a massless fluid (such as a phonon gas) but it would act on 
the thermal kaons. 

A similar problem that we need to consider is that of superfluid neutrons and protons 
coupled to thermal excitations [24, 65]. In the simplest approximation one may assume that 
the system is well below the superfluid transition temperature in which case thermal effects 
can safely be neglected. This is, indeed, the assumption in most existing work on neutron 
star dynamics. However, given the strong density dependence of the superfluid pairing gaps 
there will always be regions of the star that are close to the superfluid transition temperature. 
In these regions thermal effects are important. To model such regions one should consider a 
three-fluid system of superfluid neutrons, superconducting protons (locked to the electrons) 
and thermal excitations (the entropy). Formally, such a system would be very close to the 
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CFL-K° core. It may be sufficient to replace the density of the CFL condensate with the 
neutron density and the kaon density with the proton density in all the above equations. 

VIII. CONCLUSIONS 

We have presented a general flux- conservative formalism for modelling dissipation in 
multi- fluid systems, extending and correcting the model from [4]. The formalism was devel- 
oped with superfluid neutron star cores in mind, but is sufficiently general that it can be 
applied to a variety of analogous multi-fluid systems. The introduction of extra degrees of 
freedom (the separate "fluids" ) leads to a number of additional dissipation coefficients com- 
pared to standard single-fluid hydrodynamics. This may affect the dynamics of the system 
significantly, e.g. impact on the nature of the modes of oscillation of a neutron star. This 
may, in turn, have repercussions for the gravitational waves emitted by the system. In order 
to understand this effect we developed the most general form for the dissipation coefficients 
and discussed how their number can be reduced by imposing the "superfluid" constraint 
on one or more fluids (imposing that the flow must remain irrotational) . This constraint is 
often too severe, as a superfluid condensate will not be macroscopically irrotational, but will 
rather mimic the effect of bulk rotation by creating an array of vortices. This leads to new 
dissipative effects, such as the vortex mediated mutual friction. The superfluid constraint 
does, however, drastically simplify the dissipation problem and may be a reasonable approx- 
imation at temperatures well below the superfluid transition temperature. The simplifled 
model may also give some insight into which dissipation channels are likely to make the 
most important contributions. 

As examples of relevant applications, that take us (quite far) beyond the level that has 
been previously considered, we outlined two problems where a three-fluid description is 
appropriate. In fact, although two-fluid hydrodynamics is a useful approximation for the 
study of a system of neutrons and protons at low temperature (i.e. the outer core of a 
neutron star) [4, 58], there will always be regions of the star that are close to the superfluid 
transition temperature and in which the thermal excitations of the system (the entropy in our 
language) should be considered as a separate component [24]. Furthermore, several equations 
of state predict the appearance of hyperons in neutron star cores, with A and S~ expected 
to have the lowest threshold densities [11, 62]. The effects of multi- fluid hydrodynamics on 
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the bulk viscosity damping of the gravitational- wave driven r-mode instability have already 
been studied in [6] , although in that case the model was simphfied to the two-fluid level by 
neglecting the A hyperons. A "complete" model would require a three-fluid description to 
account for neutrons, A hyperons and a charge neutral fluid of protons, electrons and E~ 
hyperons. 

Another interesting possibility is that the ground state of matter at the extreme densities 
of neutron star cores may correspond to deconflned quarks. At asymptotically high densities 
matter is expected to be in the CFL state, and one can thus model the system by considering 
a two-fluid model of a CFL condensate coupled to phonons (the entropy) . For realistic core 
densities it is, however, believed that the ground state will be represented by the so-called 
CFL-K'' phase [21] in which one also has a kaon condensate. Motivated by this, we presented 
an example of a three-fluid system given by two superfluids, the CFL condensate and the 
kaons, coupled to a phonon gas. This example is interesting also because it can easily be 
adapted to describe neutrons and protons coupled to thermal excitations by replacing the 
kaon and CFL densities with the neutron and proton densities. The two systems are formally 
equivalent. 

From a technical point-of-view, the model developed in this paper takes us to the level 
where we need to focus on more detailed applications. In order to do this, we must give 
some thought to the nature of the many new dissipation coefficients that the system allows. 
Are there situations where the additional channels are important? Our experience from 
neutron star oscillations and the case of the superfluid mutual friction suggests that the 
answer is non-trivial (as one would expect), and that one may have to work out the detailed 
dynamics of each model system before knowing for certain. A better insight into the relative 
strength of the different dissipation coefficient and their nature (e.g. the scaling with the 
key variables of the problem, like density and temperature) would be very useful, and may 
suggest suitable simplifications for any problem at hand. Of course, this presents us with 
a challenge given that most considerations of NS transport phenomena are based on the 
single-fluid model. We need to move beyond this level to make further progress. 
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